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Аннотация 


Введение. Управление вибрационными взаимодействиями актуально в производственных процессах, связанных 
с жидкими и сыпучими средами, системами твердых тел, испытывающих кинематические или силовые 
возмущения. При этом нет единой методологической основы для формирования вибрационных взаимодействий. 
Не решены вопросы построения оптимальных вибрационных полей технических объектов. Цель 
исследования — развитие структурного подхода к разработке математических моделей в задачах формирования, 
оценки и коррекции вибрационных полей технических объектов в условиях интенсивных силовых и 
кинематических нагружений. Ставится задача построения вибрационных полей, оптимальных по совокупности 
требований, с возможностью выбора критерия оптимальности вибрационного поля технического объекта. 
Материалы и методы. В качестве базовой методологии применяется структурный подход. Он основан на 
сопоставлении механических колебательных систем, используемых как расчетные схемы технических объектов, 
и структурных схем систем автоматического управления, эквивалентных в динамическом отношении. Для 
структурного математического моделирования использовали формализм Лагранжа, элементы операционного 
исчисления на основе интегральных преобразований Лапласа, разделы теорий колебаний, алгебраические 
методы, теорию сплайн-функций. 

Результаты исследования. Предложен подход к выбору критериев оптимальности вибрационных полей на 
основе минимизации невязки вибрационных полей для различных необходимых условий. Проблема 
рассматривается в рамках механической колебательной системы, образованной твердыми телами. Показано, что 
оптимальное вибрационное поле определяется внешним возмущением и должно удовлетворять условию Ау = В. 
Здесь А — матрица, отображающая оператор условий на форму вибрационного поля в контрольных точках; 
Ь — вектор значений характеристик вибрационного поля; «-—» над у означает амплитуду колебания 
установившейся компоненты координаты. Для оценки поля с учетом зашумленных или недостоверных 
требований к динамическим характеристикам используется параметр сглаживания, обозначающий приоритет 
критерия оптимальности формы вибрационного поля. Построение поля для механической колебательной 
системы показало, что значение амплитуд колебания обобщенных координат сохраняется постоянным при 
изменении частоты внешних кинематических возмущений. Рассмотрены два подхода к коррекции критериев 
оптимальности поля: уравнивание амплитуд колебаний координат технического объекта и выбор 
энергетического оператора. 

Обсуждение и заключение. Развитие прикладной теории оптимальных вибрационных полей предполагает, во- 
первых, сопоставление оператора энергии и оператора требований к форме вибрационного поля в теории 
абстрактных сплайнов. Вторая пара сопоставляемых элементов — критерий оптимальности вибрационного поля 
и система требований к характеристикам поля в контрольных точках. Усовершенствованная таким образом 
структурная теория оптимальных вибрационных полей найдет применение в разных отраслях. Точные расчеты в 
формировании, оценке и коррекции состояний систем при вибрационном нагружении необходимы в задачах 
повышения долговечности конструкций, улучшения измерений в сложных колебательных системах, разработки 
новых технологий и материалов. 
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Введение. Особенности вибрационных взаимодействий учитываются при создании производственных 
технологий [1]. Использование вибрационных методов интенсификации технологических процессов требует 
развития математического моделирования динамических состояний технических объектов в условиях 
вибрационного нагружения. В научной и прикладной литературе прорабатываются вопросы, связанные с теорией 
колебаний и автоматического управления, теоретической механикой и динамикой машин [2]. В частности, в 
динамике машин особенности вибрационных взаимодействий учитываются в двух типах задач. Первый: задачи 
виброизоляции и виброзащиты [3]. Второй: задачи, связанные с оценкой, формированием и коррекцией 
состояний технических объектов в условиях интенсивных силовых или кинематических нагружений [4]. В 
последнем случае обращают на себя внимание эффекты вибрационных взаимодействий, которые находят 
применение в таких технологических процессах, как: 

— транспортировка горных пород; 

— разделение сыпучих материалов на фракции; 

— упрочнение поверхностей длинномерных элементов авиационных конструкций; 

— пространственная ориентация деталей в процессе сборки. 

Задачи виброзащиты связаны с решением проблем безопасности железнодорожного транспорта. Речь идет, в 
частности, о вопросах динамического состояния элементов подвески [5], шуме и вибрации [6]. Задачи оценки 
динамических состояний, вызванных значительными нагрузками, возникают при оценке вибрационных 
характеристик рельсового пути с композитными шпалами [7]. Требуют контроля вредные колебания, 
распространяющиеся в городской среде от движущегося состава, и а этом случае тоже необходима оценка 
сложных динамических состояний [8]. Аналогичные примеры связаны с отслеживанием взаимодействий 
тележки вагона с рельсами [9]. Динамическое состояние формируется при нарушении скользящего контакта 
головки пантографа с контактным проводом [10]. В данном контексте следует упомянуть и задачи анализа 
реакции рельса на вертикальные воздействия со стороны движущегося транспортного средства [11]. 

Еще одна обширная группа задач, требующих оценки динамических состояний элементов, — вибрационные 
процессы технологических машин [12]. Так, динамическое состояние машин вибрационного грохочения зависит 
от формирования динамических состояний нескольких корпусов, которые можно рассматривать как твердые тела 
с учетом характера упругих связей между ними [13]. Существенный фактор, определяющий динамические 
состояния объектов технологического назначения, — частота колебаний [14]. Отдельно необходимо упомянуть 
вибрационные транспортные машины, динамическое состояние которых определяет особенности движений [15]. 

Стоит отметить результат обобщения задач оценки, формирования и коррекции динамических состояний 
технических объектов в условиях интенсивных нагружений. Данный подход открыл возможность уточнить 
представления о том, как распределения амплитуд колебаний координат контрольных точек зависят от связности 
внешних возмущений. 

Параметры внешних вибрационных возмущений определяют возможности реализации динамических 
эффектов, которые могут действовать разнонаправленно: 

— представлять опасность для функционирования технических объектов; 

— повышать производительность технологических процессов. 

Пример негативного эффекта — превышение допустимых значений амплитуд колебаний. Полезные эффекты: 
существенное уменьшение амплитуд колебаний координат и установление определенного распределения амплитуд. 

В задачах вибрационного упрочнения длинномерных деталей под динамическими эффектами можно 
понимать совпадение форм колебаний контрольных точек рабочих органов вибрационных технологических 
машин. Это обеспечивает одномерность и однородность вибрационного поля, что согласуется с требованиями к 
динамическому качеству взаимодействий рабочей среды и поверхности упрочняемой детали. Необходимую 
форму движения рабочего органа можно выразить с помощью условий равенства амплитуд колебаний 
контрольных точек. Для исследования понятия вибрационного поля обобщается понятие динамического 
состояния в точке — одновременно рассматривается совокупность точек, распределенных по некоторому 
многообразию точек твердого тела или системы твердых тел [16]. 

В задачах вибрационных взаимодействий можно отобразить ряд динамических эффектов с помощью 
рычажных, упругих, межпарциальных неудерживающих связей. Все они — факторы специфических 
динамических эффектов, которые проявляются под воздействием внешних возмущений особенной 
конфигурации. Для определения таких конфигураций внешних возмущений, обеспечивающих динамические 
эффекты, предполагается разработка методологической базы, позволяющей строить математические модели 
вибрационных полей с учетом системы требований к особенностям взаимодействий элементов системы. 
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Предполагаемая методологическая база должна обладать: 

— эффективным математическим аппаратом; 

— инструментом наглядного отображения структурных особенностей; 

— возможностями построения динамических аналогий для обобщения разработанных методов на широкий 
класс физических процессов. 

Материалы и методы. Методологической базой может служить структурный подход, основанный на 
сопоставлении механических колебательных систем и систем автоматического управления [17]. В этом случае 
механическая колебательная система рассматривается как расчетная схема технического объекта. Ее 
сопоставляют со структурной схемой эквивалентной в динамическом отношении системы автоматического 
управления. Установленное соответствие позволяет анализировать и оценивать динамические характеристики 
технических объектов с учетом представлений об установившихся формах колебаний механической системы на 
основе аппарата теории автоматического управления [18]. У технологии структурного математического 
моделирования есть особенности, связанные с оценкой динамических состояний. Для нее используют 
передаточные функции структурных схем. Входные сигналы для названных функций — внешние возмущения, 
выходные — колебания координат оцениваемого объекта. Каждую координату системы можно соотнести с 
показателем амплитуды установившегося колебания. Это позволяет рассматривать распределение амплитуд по 
точкам системы как вибрационное поле технического объекта. Построение вибрационного поля с учетом ряда 
требований относится к классу обратных задач динамики, которые в общем случае могут быть некорректными. 
Один из методов решения некорректных задач связан с регуляризацией исходной задачи путем введения 
критериев оптимальности. 

Перечисленные особенности оценки динамических состояний технических объектов указывают на 
необходимость разработки общего подхода в методологии построения вибрационных полей с учетом критериев 
оптимальности. Вместе с тем методология формирования, оценки и коррекции вибрационных полей с заданными 
характеристиками недостаточно детализирована. В частности, не предусмотрена возможность строить 
вибрационные поля по критериям оптимальности. 

Представленная работа посвящена развитию методологии структурного математического моделирования. 
Для этого решается класс частных задач построения вибрационных полей технических объектов с условиями по 
амплитудам колебаний, которые отображают требования к динамическим эффектам и учитывают критерии 
оптимальности. 

Основные положения. Постановка задачи. Рассматривается модельная механическая колебательная 
система, образованная совокупностью твердых тел в условиях вибрационного нагружения кинематической 
природы. Модельная система (рис. 1) имеет 5 степеней свободы, состоит из 4 звеньев в виде твердотельных 
стержней с массами М, моментами инерции „Л, центрами тяжести в точке О, отстоящих от крайних 
точек А; и А: на расстояния ги р, где! = 1...4. 

Твердые тела в виде стержней сочленены между собой в крайних точках — это места крепления упругих 
элементов, установленных на опорных поверхностях, совершающих синфазные гармонические колебания (рис. 1). 


У А У? А У? Аз У А4 У5 
а О! 2 бы 0: 2 2 ы О: 5 2 г © 142 — 
(А 


ен) 
7 Л 201 М! № Л 202 М> | Л 203 М3 |7 Ил 204 М4 5 


21 22 23 24 25 


Рис. 1. Механическая колебательная система, образованная твердыми телами А! 42, А2 Аз, Аз А4, Ад А5: 
А — точки сочленения твердых тел; О; — центры тяжести; 2; — опорные поверхности; „Л — моменты инерции 
Под воздействием внешних кинематических возмущений 2; со стороны опорных поверхностей система 
совершает малые установившиеся колебания относительно положения статического равновесия. Внешние 
кинематические возмущения представляют собой связные синфазные колебания: 


Ее. (1) 
Здесь 2, — гармоническое возмущение вида 451т(0%), где © — частота и у,— коэффициенты связности внешних 
возмущений; / = 1...5. В общем случае коэффициенты связности (1) — это внешние возмущения определенной 


структуры. 
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Необходимые динамические эффекты обеспечиваются совокупностью модельных условий, а именно: 
амплитуды колебаний координат контрольных точек технического объекта принимают фиксированные значения 
вне зависимости от частоты внешних возмущений. 

Задача — определить для механической колебательной системы такие внешние кинематические возмущения, 
чтобы соответствующее вибрационное поле в контрольных точках удовлетворяло заданной совокупности 
условий в виде фиксированных амплитуд колебаний. 

Результаты исследования 

Математическая модель. Для составления системы дифференциальных уравнений рассматриваются 
координаты {ут, уз, уз, 4, Уз} смещений точки А, 1 = 1...5 относительно положений статического равновесия. 
Наравне с координатами {у1, уд, уз, ул, у5} рассматриваются: 

— координаты {201, 202, 203, 204, 205}, отображающие смещения центров тяжести О; стержней относительно 
положений статического равновесия; 

— координаты {Фу, фо, фз, фл, Ф5} малых углов поворота вокруг центров тяжести. 

Координаты 2%, Фь у; связаны между собой выражениями: 

20: = @у: + У (2) 
ф: =сКун —#/” 
[2 [а 1 


здесь а; = ‚= ‚С = . 
1 + [2 1 +72 1 +12 


Математическая модель в виде дифференциальных уравнений Лагранжа 2-го рода строится на основе 
выражений для потенциальной и кинематической энергии. Их можно представить в матричном виде с помощью 
скалярного произведения (,) векторов в А": 


1 =& = - = 1 ы - 1 = = 
П=5(К-2/5 2), Т=5 (М, )+ (19.5), (3) 
здесь 2о,ф — производные координат 20, ф ; 
К 0 Л 0 М! 0 
К У М 
= 2 = 2 в 2 
Кз Уз Мз 
0 К4 0 /4 0 М4 
У 21 
Ф1 201 
У2 22 
5 ф>2 = 202 5 = 
ф= ‚20 = ‚У=| Уз 2=| 23 
фз 203 
У4 24 
Ф4 204 
У5 25 


и = 4 = 5 5 
Векторы 20 е “, фЕЁ*, УЕК” сучетом (2) связаны между собой соотношениями: 


20 = (У, ф=У>. (4) 
Ч Ы 0 —С С 0 
Ь о 
Здесь И = — ‚Г = м. 
аз 6з —с3 сз 
0 Аа4 Ь4 0 —С4 С4 


В координатах {у}! кинетическая энергия (3) принимает вид: 
— 1 у у т 
т=5 (о МИ+У л/}5,5). (5) 
С учетом выражений (2)—(5) уравнения Лагранжа 2-го рода принимают вид: 
(ИТМИ УТ) 5+ Ку = КЕ. (6) 


Предполагается, что начальные условия — нулевые. В этом случае интегральные преобразования Лапласа 
приводят (6) к системе алгебраических уравнений: 


(Мор?+К)у= КЕ, (7) 


где Мо=ИТМИ+УГ.Л/ — матрица системы, знак «—» над векторами у и 2 обозначает преображение по Лапласу [19]. 
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С помощью методов [18] алгебраическую систему (7) можно представить в виде структурной схемы 
эквивалентной в динамическом отношении системы автоматического управления (рис. 2). 
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Рис. 2. Структурная схема механической колебательной системы (рис. 1): р=/, /=\-1 


На основе структурной схемы для фиксированного набора коэффициентов связности у; можно построить 
передаточные функции. Они определяются выбором объекта, динамическое состояние которого оценивается по 
выражениям: 


РЕ, (р) = (р) = 2 (р) =, (р) =>. (8) 
28 28 28 28 28 

Амплитудно-частотные характеристики передаточных функций (8) в физическом смысле рассмотрим как 
рычажные СВЯЗИ, представляющие собой отношения амплитуд колебаний координат Уд К амплитуде колебаний 
внешнего кинематического возмущения {2%}. 

В рамках рычажной интерпретации положительные ветви амплитудно-частотной характеристики 
отображают рычажные связи 1-го рода (реализуются двуплечими рычагами). Отрицательные ветви отображают 
рычажные СВЯЗИ 2-го рода (реализуются двуплечими рычагами). Режимы резонанса и динамического гашения 
колебаний отображают вырожденные рычажные связи, которые можно представить условными виртуальными 
рычагами [© нулевыми ИЛИ бесконечными плечами. 

Выбор структуры внешних кинематических возмущений путем задания коэффициентов связности {у 
однозначно определяет амплитуды колебаний координат системы. Вместе с тем требования, предъявленные К 
амплитудам колебаний, могут: 

== предполагать неоднозначность решения, 

— прямо противоречить друг другу; 

= допускать лишь частичное соблюдение из-за больших погрешностей. 

В таких ситуациях необходимы критерии выбора внешних возмущений, которые однозначно определят 
вибрационное поле по совокупности требований. Логичный ДЛЯ этой задачи подход основан на том, что внешние 
кинематические возмущения, обеспечивающие необходимые условия, можно определить по критерию 
оптимальности вибрационных полей. Суть данного критерия: амплитуды колебаний технического объекта 
подбираются таким образом, что доставляют минимум определенному энергетическому функционалу [20]. 
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Формирование вибрационных полей механических колебательных систем на основе критерия 
оптимальности. В рамках модельной задачи предполагается, что в точках Ая, Ар...А вибрационного поля {уд, 
1=1...п механической колебательной системы задано определенное условие. Оно заключается в том, что 
амплитуды {уж}, А = 1...т принимают фиксированные значения {уо1...Уот!: 

Ук = Уи. (9) 

Здесь А = 1...т, где т отображает количество точек, в которых заданы условия на амплитуды. В этом случае 
и далее по тексту знак «—» наду или 2 означает амплитуду колебания установившейся компоненты координаты, 
если не оговорено иное. В остальных точках {у» {= 1...п}/{ун, К = 1...тф не заданы требования к значениям 
вибрационного поля. 

Необходимые требования (9) можно представить в виде: 

АУ=Ь, (10) 


где А — матрица, отображающая оператор условий на форму вибрационного поля в контрольных точках; 
р — вектор значений характеристик вибрационного поля [20]. 

В общем случае условия для вибрационного поля (10) могут быть неопределенными или несовместными. Для 
рационального учета несовместных или неопределенных требований к вибрационному полю вводится критерий его 
оптимальности и допущение, что совокупность требований может быть выполнена с некоторым приближением. 

Под значениями у1, у2...у„ понимаются амплитуды колебаний установившихся компонент, поэтому для 


каждой фиксированной частоты ® внешних возмущений 2 можно поставить вариационную задачу построения 
оптимального вибрационного поля, которое: 

— доставляет минимум энергетическому функционалу; 

— отображает представления о формах колебания или о значениях потенциальной или кинетической энергии. 

В общем случае энергетический функционал определяет особенности распределения амплитуд динамических 
характеристик по точкам механической колебательной системы. 

В рамках рассматриваемой модельной задачи энергетический функционал, заданный оператором Те, 
отображает неравномерность вибрационного поля. Ее показатель — сумма квадратов разностей отклонений 
амплитуд колебаний соседних точек: 


п-1 


т.(у)= У (у -Уь. (0 
1=1 


Нулевое значение энергетического функционала (11) определяет равномерное вибрационное поле, амплитуды 
колебаний точек которого равны между собой. Минимизация энергетического функционала с учетом заданных 
условий на амплитуды колебаний в контрольных точках (10) определяет форму вибрационного поля, 
экстремально близкую к равномерной. 

Примем во внимание ряд условий: 

— фиксированная частота внешних возмущений 0; 

— вибрационное поле оптимально в смысле критерия, заданного энергетическим функционалом; 

— амплитуды вибрационного поля в контрольных точках удовлетворяют условиям (10). 

В этом случае формальную постановку задачи построения вибрационного поля можно представить в виде: 


Ув(®)= аге тт Т.( У), (12) 


Ау=Ь,у=(-Моф?+К) К 
где у — вибрационное поле, заданное внешними кинематическими возмущениями 2 , которые рассматриваются 


как независимые вариационные параметры. 

Для фиксированной частоты © вибрационное поле уо(0), точно удовлетворяющее совокупности требований, 
формально находится как решение (12) условной минимизации квадратичной функции (11) с учетом линейных 
условий (10): 


Цу, Л) = (Тоу,у)+(АУ-Ь, Л) > тт |. (13) 
л,у=(Мор?+К) К 
Здесь у — вибрационное поле, определяемое внешним кинематическим возмущением; 
Л = (41...\.) — коэффициенты Лагранжа. Минимизация в задаче (13) идет по параметрам = внешних 


возмущений. Для фиксированных значений 2 амплитуд внешних возмущений создаваемое вибрационное 
поле У является решением задачи: 


Мо? +К)у = КЕ, (14) 


где Мо=И"ТМИ +ИГЛ/ — матрица системы. 
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Вибрационное поле у, точно удовлетворяющее заданным значениям (10), можно определить на основе 


решения задачи условного минимума с помощью функции Лагранжа: 
—1 —1 —1 
ЦЕ, Л)= (т (Мор?+К) КЕ,(Мор?+К) к) (4(мор? +К) КЕ-Ьл) —> тт, (15) 
Ат 


где = — вектор амплитуд колебаний внешних кинематических возмущений; То — матрица, определяющая 
квадратичную форму энергетического функционала (11). 
Внешние возмущения 2., формирующие необходимое вибрационное поле, множители Лагранжа Л 
находятся решением системы: 
т 2 иг 2 = 2 Е" 
КТ |(-Моо?+К) | То(-Мов?+К) К [4(-Моо?+К) К] |(=)\ (0 
=| |. (16) 
ь — Л Ь 
А(-Мос?+К) К 0 
Итак, 2‹ — это решение задачи условной минимизации (15), (16). Оптимальное вибрационное поле ус 
определяется внешним возмущением 2. и имеет вид: 


уз = (-Моф?+К)" Кго. (17) 


Здесь ус можно назвать интерполирующим вибрационным полем, так как при условии совместности (16), оно 
точно удовлетворяет требованиям (10). С учетом зависимости от частоты , совместность и определенность 
системы (16) характеризуется в рамках теорем существования и единственности решений систем уравнений 
теории сплайнов [20]. При условии совместности и определенности системы (16) и (17), оптимальное 
вибрационное поле у. точно удовлетворяет совокупности условий (10). 

Необходимые условия (10) могут содержать зашумления, ошибки или быть противоречивыми. С учетом 
противоречивых или недостоверных требований, вибрационное поле можно построить с помощью решения 
задачи сглаживания. Ее результат — оптимальное вибрационное поле, удовлетворяющее условиям, отличным от 
исходных требований, но близких к ним. 

Оценка вибрационного поля с учетом зашумленных или недостоверных требований к динамическим 
характеристикам. Вибрационное поле у(®) (17), удовлетворяющее условиям (10), может быть 
аппроксимировано с помощью вибрационного поля уо(6) с учетом недостоверности или противоречивости 
необходимых требований (10): 


у(Ф) = ты, | 1 у)+У и -уи |. (18) 


7=(-Мов?+к) К = 
Здесь Теу — оператор, отображающий критерий оптимальности формы вибрационного поля; а — параметр 


сглаживания, обозначающий приоритет критерия оптимальности формы вибрационного поля по отношению к 
требованиям выполнения условий (10). 

В выражении (18) 2 рассматривается как независимая переменная. Вибрационное поле у.(©) можно назвать 
сглаживающим, так как оно приближенно удовлетворяет системе необходимых требований. 

Вибрационное поле у, определяется системой уравнений, полученной путем минимизации функции 
Лагранжа: 


1. (У(Е)) = (Тау, )+М-Ы[ > тт. (19) 
Здесь У(2) является вибрационным полем, формируемым внешним кинематическим возмущением 27: 
(-Моо? +К)у = КЕ. (20) 
Сделаем замену на основе (20): 
У= (-Мов? +К) КЕ. (21) 


Теперь функция Лагранжа (19) с учетом (20) и (21) приводится к виду: 
= 2 


№ (7 )= сто (-моо?+К) ' К2.(-Мов? +К) КЕ)+ 4(-Мнь: +К) КЕ-Ь 


> тт, (22) 


где | — \/(.,-) —_ длина вектора в В", выраженная через скалярное произведение (. .) . 
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Рассмотрим внешние кинематические возмущения 2. Они формируют оптимальное вибрационное поле у, 
которое сглаживает заданные условия (10) и является решением задачи минимизации (22). Такие возмущения 
определяются системой: 


Т Т 
= = =1 
т (Мо? +К) (яТо + АТА) (-Мов? +К) Кла =К' (Мо? +К) АТЬ. (23) 
После упрощений системы (23) искомые амплитуды кинематических возмущений 2. определяются 
выражением: 
= 2 ва тат 

2а =К '(Мор?+К)(оТь+ АТА) АТЬ. (24) 
Как показано выше, внешние возмущения 2. (24) — это решение задачи (23) с учетом выбора параметра 
сглаживания ц. Эти возмущения обеспечивают построение вибрационного ПОЛЯ ус, кОоТОрое приблизительно 

удовлетворяет совокупности требований (10) с учетом возможной несовместности или недостоверности: 


2 —1 
У« =(Мор?+К) Кга. (25) 
Вид вибрационного поля ух (25) явно не зависит от частот внешних возмущений ©: 
—1 
Уа = (а7То+ АТА) АТ. (26) 


Внешнее возмущение 2, обеспечивающее построение оптимального сглаживающего вибрационного поля 
Уз (24), можно использовать для построения интерполяционного вибрационного поля ус по мере «0. Выбор 
параметра сглаживания @ определяется из условий малости невязки: 


2 
ф(о,) = [Ау -В, (27) 
где у = (-Моо? + к)" К2а, 2а — кинематическое возмущение (27). 


Необходимо отметить, что вибрационное поле у» зависит от частоты © неявно, через параметр сглаживания 
а = (6), выбор которого определяется частотой внешнего возмущения 6 и значением невязки (27). 

Построение вибрационного поля для механической колебательной системы. Вычислительный 
эксперимент. Пусть параметры механической колебательной системы (рис. 1) заданы жесткостями К», 
массами Мь, моментами инерции Л. и геометрическими характеристиками Пь, 5; (таблица 1). 


Таблица 1 
Параметры модельной механической колебательной системы 
№ М, кг У, кг-м? К, Н/м 1, м Ь, м 
1 100 10 1000 1 1 
2 200 20 2000 1 1 
3 300 30 3000 1 1 
4 400 40 4000 1 1 


В качестве требований к амплитудам вибрационного поля заданы условия на значения в точках А; (таблица 2). 


Таблица 2 
Амплитуды колебаний в точках вибрационного поля (у, м) 
Значения А 42 Аз Ал А5 
Необходимые 0,002 — — -- 0,005 
Сглаженные 0,002 0,0027 0,0035 0,0043 0,005 


В качестве энергетического функционала используется выражение (11), которое можно представить с 
помощью оператора То: 


Я ЕЕ я 
—_ у2- Я 1 —1 0 —_ 
У? ЕЕ У? 
— Уз—у2 1-1 — 
Да Уз = = |= Уз 
_ 74 — Уз 1-1 _ 
У4 = У4 
— Уз У4 0 Е 
У5 У5 
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Выражение (7. су, 1 .У) можно представить в виде: 


= о\( я) (у 
12 - У || у 
(тьу,Тьу) = (ть Тьу, у) = (Ту, у) = 12 -1 Уз || Уз 
1 2-1 94 || У4 
0 ГОТ 5) (Уз 
Здесь То — матрица, определяющая квадратичную форму. Значение оператора Ту равно нулю на 


равномерном вибрационном поле у = соп5!. Амплитуды колебаний точек этого поля равны между собой для 


внешнего кинематического возмущения на произвольной частоте. 
Зададим условия (10) на вибрационное поле. Для этого можно использовать оператор данных А15, который 


сопоставляет вибрационному полю у значения амплитуд в точках Ат, 45: 


я я 
0000] (я 
| бо, 

У У 

У5 У 


Вибрационное поле, интерполирующее заданные значения в точках А, 45 (таблица 2), можно 
аппроксимировать с помощью сглаживающего вибрационного поля у» для параметра сглаживания а = 0,01. 
Найдем разницу между значениями вибрационного поля и требуемыми значениями для выбранного 
коэффициента сглаживания @. Она определяется функцией невязки (27), значение которой составляет 
приблизительно 0,00001 м, где Б= (0,002; 0,005). Построенное вибрационное поле, сглаживающее 
интерполяционные условия, имеет линейную форму (таблица 2). Для фиксированного параметра сглаживания 
а = 0,01 совокупность амплитуд колебаний точек вибрационного поля у. обеспечивающего близость к 
интерполируемым данным с невязкой ф(а) = 0,00001 м (таблица 2), имеет вид постоянных функций частоты 


внешнего возмущения в (рис. 3 а). 


(4) (5) 2, м (5) 


10 о, рад/с 
(0) 


6) 
мы (5) Ф(©), м 
4 
© 0,0003 
0,006 (3) 
0,002 (1) (2) 0,0001 
0 1 2 3 4 5 
п, номер обобщенной координаты 0 2 4 6 ®, рад/с 
в) г) 


Рис. 3. Характеристики вибрационного поля: а — значения амплитуд колебаний сглаживающего вибрационного поля для 
частот внешних возмущений на интервале в Е (0,100) рад/с; 6 — амплитуды внешних кинематических возмущений, 
обеспечивающие требуемые амплитуды колебаний в зависимости от частоты возмущений 0; в — форма вибрационного 
поля для частот на интервале в Е (0...100); (1), (5) — требуемые значения вибрационного поля; (2), (3), (4) — значения, 
построенные на основе энергетического подхода; г — отклонение ф (©) вибрационного поля от необходимых значений 
амплитуд в фиксированных точках технического объекта для постоянного коэффициента сглаживания 41 = 0,1 
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Значение амплитуд колебаний обобщенных координат в данном случае сохраняется постоянным при 
изменении частоты внешних кинематических возмущений 2. Эти возмущения, необходимые для обеспечения 
требуемых амплитуд в фиксированных точках вибрационного поля, существенно зависят от частоты © внешних 
возмущений (рис. 3 0). 

Опишем вибрационное поле для фиксированной частоты или частотного интервала, на котором оно остается 
неизменным. Его можно представить графиком функции, зависящей от номера обобщенной координаты или 
координаты точки системы, в которой оценивается вибрационное поле. Форма вибрационного поля определяется 
критерием оптимальности. В частности, энергетический функционал То, определяющий критерий 
оптимальности, формирует линейную форму вибрационного поля (рис. 3 в). Необходимо отметить, что 
«интерполяционное» вибрационное поле строится с помощью «сглаживающего» вибрационного поля, 
соответствующего параметру сглаживания а. Отклонение «сглаживающего» вибрационного поля должно быть 
незначительно. В рассматриваемом модельном примере для а=0,1 величина отклонения сглаживающего 
вибрационного поля представляется константой ф (9) = 0,0001 м (рис. 3 г). Необходимо отметить, что в общем 
случае энергетического функционала значение коэффициента сглаживания а, обеспечивающего фиксированный 
уровень невязки ф(а) = Ф’, зависит от частоты @ и системы требований. 

Таким образом, варианты формирования динамических состояний механических колебательных систем в 
виде вибрационного поля, удовлетворяющего системе условий, зависят от: 

— энергетического функционала, отражающего особенности вибрационного поля; 

— оператора, отображающего совокупность требований к вибрационному полю. 

Возможны и другие варианты формирования вибрационного поля на основе выбора энергетического оператора. 

Коррекция критериев оптимальности формообразования вибрационных полей с помощью весовых 
коэффициентов энергетических операторов. Выше рассмотрена форма вибрационного поля с возможностью 
уравнивания амплитуд колебаний координат технического объекта. Есть и другие подходы к минимизации 
амплитуд обобщенных координат, например, с помощью выбора энергетического оператора Те в виде 
тождественного отображения: 

ТЕ: у> у, 
где Е — тождественный оператор. 

Формирование вибрационного поля на основе энергетического оператора Те связано с минимизацией или 
обнулением квадратичной функции амплитуд синфазных гармонических колебаний фиксированной частоты. 
В результате решения задачи сглаживания для энергетического функционала, соответствующего единичному 
оператору Ть, формируется вибрационное поле. Амплитуды его колебания в точках 4›, Аз, Ад равны 


нулю (рис. 4 а), т. е. в точках 42, Аз и А4 реализуется режим динамического гашения колебаний. 
Ко 


—0,005 


| О ПО 3 5 


] п, номер обобщенной координаты —0,015 
а) 6) 
г (5) Ум 
0,005 Мид 
0,006 
0,003 (0 
0,002 
0,001 о 9 @& ь 
0 5 4 р: ®, рад п, номер обобщенной координаты 


в) г) 


Рис. 4. Возможности коррекции вибрационного поля: а— форма вибрационного поля для энергетического функционала 
отображает минимально возможные амплитуды колебаний с учетом необходимых требований; 6 — амплитуды внешних 
кинематических возмущений 21, обеспечивающих необходимые требования к динамическим особенностям вибрационного 
поля; в — амплитуды колебаний вибрационного поля технического объекта в зависимости от частоты внешнего 
возмущения; г — варианты вибрационных полей в зависимости от весового коэффициента энергетического оператора Тв 
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Совокупность внешних возмущений кинематического характера — это решение задачи минимизации 
энергетического функционала, в общем случае заданного на динамических состояниях механических 
колебательных систем с учетом требований к амплитудам колебаний в фиксированных точках вибрационного 
поля (рис. 46). 

Особенность вынужденных кинематических возмущений, обеспечивающих требуемую форму 
вибрационного поля, — возрастание амплитуды колебаний по мере роста частоты возмущений (рис. 4 в). 
Внешние кинематические возмущения 2 обеспечивают на различных частотах ® фиксированное вибрационное 
поле технического объекта (рис. 4 г). Параметр сглаживания а = 0,01 обеспечивает величину абсолютного 
отклонения от требуемых значений на уровне ф(а) = 0,000053. Это показывает, что совокупность амплитуд 
колебаний в остальных точках вибрационного поля определяет особенности энергетического оператора на 
основе необходимых значений амплитуд вибрационного поля в фиксированных точках механической 
колебательной системы. В частности, оператор Т. направлен на формирование равномерного вибрационного 
поля, а оператор Те — на обнуление вибрационного поля. 

Можно использовать весовые коэффициенты, отображающие комбинацию критериев оптимальности. Это 
позволит совместить: 

— требования по сближению амплитуд колебаний точек вибрационных полей; 

— требования, направленные на обнуление амплитуд колебаний точек вибрационного поля. 

Варьирование весовых коэффициентов, определяющих комбинацию оптимальных критериев формирования 
вибрационного поля, фактически является способом коррекции вибрационного поля. 

Построим семейство форм вибрационных полей технического объекта, удовлетворяющих определенным 
условиям. С этой целью используем комбинированный критерий оптимальности, реализуемый с помощью 
взвешенной суммы энергетических функционалов с учетом весового коэффициента: 

|" = +В[ЬУ|. (28) 

Здесь В — весовой коэффициент, определяющий комбинацию критерия оптимальности на основе малости 
амплитуд и критерия на основе однородности амплитуд вибрационных полей. Варьирование весового 
коэффициента В приводит к построению семейства вибрационных полей. Критерий их оптимальности имеет 
сложную структуру, учитывающую разницы и абсолютные значения амплитуд колебаний точек вибрационного 
поля. Выбор совокупности весовых коэффициентов определяет соответствующие формы вибрационных полей, 
которые реализуются под воздействием внешних возмущений, являющихся решением вариационной задачи [20]. 

На основе совокупности весовых коэффициентов В можно построить оптимальные вибрационные поля, с 
достаточной точностью удовлетворяющие требованиям к амплитудам колебаний в фиксированных точках (таблица 3). 


Формы оптимальных вибрационных полей для весовых коэффициентов энергетического ао, 
№ В м Невязка Частота Е р и о а 
Ф(а), м ©, рад/с 
1 1 0,01 | 0,0000841 10 0,0019704 | 0,0009850 | 0,0009845 | 0,0019686 | 0,0049213 
2 0,5 0,01 | 0,0000503 10 0,0019847 0,0014467 | 0,0016322 | 0,0026337 | 0,0049521 
3| 0,25 0,01 | 0,0000305 10 0,0019940 | 0,0018969 | 0,0022740 | 0,0032196 | 0,0049701 
4 | 0,125 | 0,001 | 0,0000195 10 0,0020000 0,0022503 | 0,0027818 | 0,0036611 | 0,0049980 
5 | 0,0625 | 0,001 | 0,0000014 10 0,0020004 | 0,0024754 | 0,0031053 | 0,0039292 | 0,0049986 
Комбинирование оптимальных критериев формирования вибрационных полей — это способ коррекции 


динамических состояний технического объекта (рис. 4 2). Отметим, что найденные кинематические 
возмущения 2:, рассматриваемые как функции частоты ©, имеют по одному нулю (рис. 5 е). 
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Рис. 5. Характеристики оптимального вибрационного поля для весового коэффициента В = 0,25 и параметра сглаживания 
а = 0,01: а — отношение амплитуд колебаний точек оптимального вибрационного поля (уа)1 к амплитуде кинематических 
возмущений (70)1; 6 — отношение амплитуд (уд)? к (2а)2; в — отношение амплитуд (уа)з к (2а)з; г — отношение амплитуд (уа)4 
к (2а)4; д — отношение амплитуд (уа)5 к (2а)5; е — значение амплитуд внешних возмущений (га) в зависимости от частоты © 


Отношения амплитуд у/2,, для которых у: = ук) и 2: = 2К®) найдены на основе решения задачи определения 
оптимальных вибрационных полей, имеют по одному разрыву 2-го рода (рис. 5 а-9). 

Следует обеспечить заданные амплитуды колебаний фиксированных точек вибрационного поля 
механической колебательной системы. Для этого можно выбрать комбинированные энергетические 
функционалы с учетом весовых коэффициентов и, основываясь на данных показателях, определить необходимые 
внешние кинематические воздействия с учетом вариационных принципов. 

Если требования к амплитудам колебаний вибрационного поля сильно зашумлены или несовместны, можно 
построить вибрационное поле, которое с определенной точностью удовлетворяет требованиям и является 
оптимальным для комбинированного критерия, но для других условий. 

Определение оптимального критерия формообразования вибрационных полей на основе функции 
невязки. В общем случае требования к вибрационному полю могут представлять собой сложный комплекс 
совместных или несовместных условий, отражающих ограничения на амплитуды смещений, скоростей, 
динамических реакций, включая условия неудерживающего характера. 

С учетом критерия оптимальности вибрационного поля механической колебательной системы совокупность 
совместных условий допускает единственное решение в форме определенной структуры внешних возмущений. 
В случае несовместности условий на вибрационное поле критерий оптимальности фактически выступает 
методом регуляризации некорректной обратной задачи. 

Энергетический функционал определяет специфику формообразования вибрационных полей на основе 
значений амплитуд колебаний смещений, скоростей, силовых возмущений или реакций, возникающих между 
элементами системы. Энергетические функционалы могут передавать представления о максимальной 
потенциальной или кинетической энергии или служить их комбинацией с учетом весовых коэффициентов, 
выполняющих роль методов коррекции вибрационных полей. 

К условиям по значениям амплитуд вибрационных полей возможно добавить условия по характеру самих 
внешних кинематических или силовых возмущений. Они выражаются в том, что амплитуды внешних 
возмущений принимают фиксированные абсолютные или относительные значения. Необходимость задания 
требований к входным возмущениям связана с тем, что для обеспечения заданных форм вибрационных полей 
амплитуды внешних возмущений могут неограниченно возрастать по мере роста частоты. 
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Итак, представим комплексную задачу выбора: 

— энергетического функционала; 

— оператора требований к вибрационному полю; 

— дополнительного оператора требований к внешним силовым или кинематическим возмущениям. 

В рамках развития структурной теории оптимального вибрационного поля можно поставить обобщенную 
задачу его построения: 


Ус (6) = а’ тт - Те (5), 
4у=Ь,У= (Мое? К) КЕ,СЕ=Х 
где С. = Х — дополнительная связь, ограничивающая размерность внешних кинематических возмущений. 
Методология выбора эффективного оператора может базироваться на минимизации функции дополнительной 
невязки, зависящей от параметров энергетического оператора: 


Ве = ат тг | ВУв -<| ь 


где ув — решение задачи определения вибрационного поля с комбинированным критерием 
оптимальности Тв ((28) и рис. 4) и условиями А (10); В — оператор определения значений в совокупности 
дополнительных точек; с — вектор значений в дополнительных точках. 

Эту задачу можно рассмотреть как задачу выбора «естественного» критерия оптимальности [21]. Таким 
образом, представленная совокупность задач и методов решений может быть фундаментальной основой 
структурной теории построения вибрационных полей. 

Обсуждение и заключение. Разработаны элементы структурной теории оптимальных вибрационных полей 
механических колебательных систем. Основные компоненты указанной теории — это методы построения 
вибрационных полей механических колебательных систем на основе вариационных принципов. Основное 
понятие — вибрационное поле механической колебательной системы. Оно представляет собой распределение 
амплитуд колебаний обобщенных координат и других динамических характеристик механической 
колебательной системы по точкам объекта с оцениваемым динамическим состоянием. 

Предполагается, что в ряде точек механической колебательной системы задана система условий. Она 
выражается через абсолютные или относительные значения смещений, скоростей и реакций. Условия — это 
требования к характеристикам вибрационного поля механической колебательной системы. Существует также 
система условий, связанных с особенностями внешних силовых или кинематических возмущений, и она тоже 
отражает требования к совокупности абсолютных и относительных значений амплитуд кинематических или 
силовых возмущений. 

Для развития прикладной теории оптимальных вибрационных полей в рамках структурных представлений 
можно установить аналогию между двумя парами. Первая образована из оператора энергии Т и оператора 
условий А в теории абстрактных сплайнов. Вторая образована критерием оптимальности вибрационного поля и 
системой требований к характеристикам вибрационного поля в контрольных точках. 

В рамках теории вариационной сплайн-аппроксимации различают интерполяционные, сглаживающие и 
смешанные сплайны. В рамках представлений о вибрационных полях механических колебательных систем на 
основе структурных представлений роль сплайна выполняет совокупность амплитуд внешних силовых и 
кинематических возмущений, а роль интерполяционных условий выполняет совокупность требований к 
амплитудам входных или выходных сигналов. 

Можно полагать, что существуют различия между вибрационными полями технических объектов, 
находящихся в условиях внешних возмущений силовой или кинематической природы, которые рассматриваются 
как ключевой фактор формирования динамического состояния. Вместе с тем, при построении вибрационного 
поля считается естественной оценка энергетических соотношений, зависящих от обобщенных координат. Их 
значения определяются по внешним силовым или кинематическим возмущениям. Фактически нужно найти два 
вектора — обобщенных координат и внешних возмущений. 

Учтем вариационный подход, а также совокупность требований к вибрационному полю и внешним 
возмущениям. Представим задачу определения внешних возмущений, формирующих совокупность 
динамических состояний, в виде вибрационного поля с учетом системы условий. В такой интерпретации речь 
идет о задаче оптимизации в разных постановках в рамках вариационных принципов. В частности, возможны 
постановки задач, аналогичных задачам построения сглаживающего или смешанного сплайна. Это зависит от 
совместности условий на входные и выходные сигналы. Отметим существенную особенность задач построения 
вибрационных полей: частота внешнего возмущения — это параметр задачи оптимизации, а в случае с системами 
с сосредоточенными параметрами векторы амплитуд служат элементами и-мерного Декартова пространства. 
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Аналогии между задачами аппроксимации и построения вибрационных полей позволяют развивать 
методологии формирования, оценки и коррекции динамических состояний технических объектов. 

Усовершенствованная структурная теория оптимальных вибрационных полей будет востребована в 
авиационной и космической промышленности. Ее аппарат можно применить для обеспечения безопасности 
полетов и долговечности конструкций. Кроме того, данный подход логично задействовать в машиностроении, 
электронике, при разработке медицинского оборудования. Теория позволит решать задачи формирования, 
оценки и коррекции состояний технических систем и устройств в условиях вибрационного нагружения. Еще одно 
направление ее практического применения — создание и улучшение измерительных способов и средств сложных 
колебательных систем. Полагаем также, что математические модели оптимальных вибрационных полей 
представляют интерес для разработки новых технологий и материалов, в частности, материалов с высокой 
амортизацией, которые могут снизить уровень вибраций и шума в технических системах. 

В качестве основного ограничения методики необходимо указать, что рассматриваются колебания координат 
объекта, которые представляют собой установившиеся формы гармонических колебаний или колебаний 
фиксированной амплитуды. 
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